function m2 = lpcvd_h2m2(h,dx)
assert(size(h,1)==size(h,2));
L = size(h,1);
if nargin == 1
    dx = pi/L;
end
dh = [diff(h,1,1);h(1,:)-h(end,:)];
m2 = sum(sum(dh.^2))/(L^2*dx^2);
